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Abstract. The temporal communication patterns of human individuals are known to 

be inhomogeneous or bursty, which is reflected as the heavy tail behavior in the inter- 
event time distribution. As the cause of such bursty behavior two main mechanisms 
have been suggested: a) Inhomogeneities due to the circadian and weekly activity 
patterns and b) inhomogeneities rooted in human task execution behavior. Here we 
investigate the roles of these mechanisms by developing and then applying systematic 
de-seasoning methods to remove the circadian and weekly patterns from the time-series 
of mobile phone communication events of individuals. We find that the heavy tails in 
the inter-event time distributions remain robustly with respect to this procedure, which 
clearly indicates that the human task execution based mechanism is a possible cause 
for the remaining burstiness in temporal mobile phone communication patterns. 
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1. Introduction 

Recently modem information-communication-technology (ICT) has opened us access 
to large amounts of stored digital data on human communication, which in turn has 
enabled us to have unprecedented insights into the patterns of human behavior and 
social interaction. For example we can now study the structure and dynamics of large- 
scale human communication networks [H [21 |3l H] and the laws of mobility O EJ [7], 
as well as the motifs of individual behavior P El IDl CB Ca [131 E]. One of the 
robust findings of these studies is that human activity over a variety of communication 
channels is inhomogeneous, such that high activity bursts of rapidly occurring events 
are separated by long periods of inactivity [151 HSl [El [IH [El EO]- This feature is 
usually characterized by the distribution of inter-event times r, defined as time intervals 
between, e.g., consecutive e-mails sent by a single user. This distribution has been found 
to have a heavy tail and show a power-law decay as -P(t) ~ r^^ [8]. 

In human behavior obvious causes of inhomogeneity are the circadian and other 
longer cycles of our lives as results of natural and societal factors. Malmgren et 
al. [21 Hn] suggested that an approximate power-law scaling found in the inter-event 
time distribution of human correspondence activity is a consequence of circadian and 
weekly cycles affecting us all, such that the large inter-event times are attributed 
to nighttime and weekend inactivity. As an explanation they proposed a cascading 
inhomogeneous Poisson process, which is a combination of two Poisson processes with 
different time scales. One of them is characterized by the time-dependent event rate 
representing the circadian and weekly activity patterns, while the other corresponds 
to the cascading bursty behavior with a shorter time scale. Their model was able to 
reproduce an apparent power-law behavior in the inter-event time distribution of email 
and postal mail correspondence. In addition they calculated the Fano and Allan factors 
to indicate the existence of some correlations for the email data as well as for their 
model of inhomogeneous Poisson process, with quite good comparison [T2] . 

However, the question remains whether in addition to the circadian and weekly 
cycle driven inhomogeneities there are also other correlations due to human task 
execution that contribute to the inhomogeneities observed in communication patterns, 
as suggested, e.g., by the queuing models [El [2I]- There is evidence for this by Goh and 
Barabasi [22], who introduced a measure that indicates the communication patterns to 
have correlations. Recently, Wu et al. have studied the modified version of the queuing 
process proposed in [8] by introducing a Poisson process as the initiator of localized 
bursty activity [23] • This was aimed at explaining the observation that the inter-event 
time distributions in Short Message (SM) correspondence follow a bimodal combination 
of power-law and Poisson distributions. The power-law (Poisson) behavior was found 
dominant for r < tq (r > tq). Since the event rates extracted from the empirical data 
have the time scales larger than tq (also measured empirically), a bimodal distribution 
was successfully obtained. However, in their work the effects of circadian and weekly 
activity patterns were not considered, thus needing to be investigated in detail. 
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As the circadian and weekly cycles affect human communication patterns in quite 
obvious ways, taking place mostly during the daytime and differently during the 
weekends, our aim in this paper is to remove or de-season from the data the temporal 
inhomogeneities driven by these cycles. Then the study of the remaining de-seasoned 
data would enable us to get insight to the existence of other human activity driven 
correlations. This is important for two reasons. First, communication patterns tell about 
the nature of human behavior. Second, in devising models of human communication 
behavior the different origins of inhomogeneities should be properly taken into account; 
is it enough to describe the communication pattern by an inhomogeneous Poissonian 
process or do we need a model to reflect correlations in other human activities, such as 
those due to task execution? 

In this paper, we provide a systematic method to de-season the circadian and weekly 
patterns from the mobile phone communication data. Firstly, we extract the circadian 
and weekly patterns from the time-stamped communication records and secondly, these 
patterns are removed by rescaling the timings of the communication events, i.e. phone 
calls and SMs. The rescaling is performed such that the time is dilated or contracted 
at times of high or low event activity, respectively. Finally, we obtain the inter-event 
time distributions by using the rescaled timings and comparing them with the original 
distributions to check how the heavy tail and burstiness behavior are affected. As the 
main results we find that the de-seasoned data still shows heavy tail inter-event time 
distributions with power-law scalings thus indicating that human task execution is a 
possible cause of remaining burstiness in mobile phone communication. 

This paper is organized as follows. In Section [2} we introduce the methods for de- 
seasoning the circadian and weekly patterns systematically in various ways. By applying 
these methods the values of burstiness of inter-event time distributions are obtained and 
subsequently discussed. Finally, we summarize the results in Section [3] 

2. De-seasoning analysis 

We investigate the effect of circadian and weekly cycles on the heavy-tailed inter-event 
time distribution and burstiness in human activity by using the mobile phone call (MFC) 
dataset from a European operator (national market share ~ 20%) with time-stamped 
records over a period of 119 days starting from January 2, 2007. The data of January 
1, 2007 are not considered due to its rather unusual pattern of human communication. 
We have only retained links with bidirectional interaction, yielding = 5.2 x 10^ users, 
L = 10.6 X 10^ links, and C = 322 x 10^ events (calls). For the analysis of Short Message 
(SM) dataset, see Appendix. 

We perform the de-seasoning analysis by defining first the observable. For an 
individual service user i, ni{t) denotes the number of events at time t, where t ranges 
from seconds, i.e. the start of January 2, 2007 at midnight, to Tf ^ 1.03 x 10^ 
seconds (119 days). The total number of events Si = X]t=o"'i(^) called the strength 
of user i. In general, for a set of users A, the number of events at time t is denoted by 



Circadian pattern and burstiness 



4 



nA{t) = J2ieA^ii't)- ^ can represent one user, a set of users, or the whole population. 
When the period of cycle T is given, the event rate PA,r(^) with < t < T is defined as 

PA,T(t) = — Mt + kT), SA = Y,Mt)- (1) 

k=0 i=0 

For convenience, we redefine the periodic event rate with period T as pa,t(^) = 
Pa,t(^ + kT) with any non- negative integer k for < t < oo. By means of the event 
rate, we define the rescaled time t*(t) as following [I2] 

t*{t) = J2 /^A,r(t')- (2) 

0<t'<t 

This rescaling corresponds to the transformation of the time variable by p*{t*)dt* = 
p{t)dt with p*{t*) = 1. Here p*{t*) = 1 means that there exists no cyclic pattern in the 
frame of rescaled time. The time is dilated (contracted) at the moment of high (low) 
activity. 

In order to check whether the rescaling affects the inter-event time distributions and 
whether still some burstiness exists we reformulate the inter-event time distributions 
by using rescaled event times and compare them with the original distributions. The 
definition of the rescaled inter-event time from the rescaled time is straightforward. 
Considering two consecutive events of a user i G A occurring at times tj and t^+i, the 
original inter-event time is r = t^+i — tj, then the corresponding rescaled inter-event 
time is defined as follows 

r*^t*it,^,)-t*{t,)= P^^Tit'). (3) 

To find out how much the de-seasoning affects burstiness, we measure the burstiness of 
events, as proposed in [22], where the burstiness parameter B is defined as 

B ^ (4) 

ar + rrir 

Here cTt and m,- are the standard deviation and the mean of the inter-event time 
distribution -P(t), respectively. The value of B is bounded within the range of [—1, 1] 
such that i? = 1 for the most bursty behavior, 5 = for neutral or homogeneous Poisson 
behavior, and B = —1 for completely regular behavior. The burstiness of the original 
inter-event time distribution, denoted by Bq, is to be compared with that of Bt of the 
rescaled inter-event time distribution for given period T. With the de-seasoning the 
burstiness is expected to decrease and here we are most interested in by what amount 
the burstiness decreases when using T = 1 day or 7 days, i.e. removing the circadian or 
weekly patterns. 

2.1. Individual de-seasoning 



First, we perform the de-seasoning analysis for individual users with various values of 
T. For some sample individuals in the MFC dataset, we obtain the original and the 
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Figure 1. De-seasoning MPC patterns of individual users: the original and the 
rescaled event rates with period of T = 1 day (left) and the original and the rescaled 
inter-event time distributions with various periods of T (right). Individual users with 
the strength Si = 200 (a), 400 (b), 800 (c), 1600 (d), and 3197 (e) are analyzed. The 
original inter-event time distribution of the whole population is also plotted as a dashed 
curve for comparison. 



rescaled event rates. The event rates in case of T = 1 day are depicted in the left 
column of Fig. [TJ The strengths of individual users are 200, 400, 800, 1600, and 3197. 
In most cases we find the characteristic circadian pattern, i.e. inactive nighttime and 
active daytime with one peak in the afternoon and another peak in the evening. The 
rescaled event rate successfully shows the expected de-seasoning effect, i.e. p*{t*) = 1, 
except for weak fluctuations. 

The rescaled inter-event time distributions up to T = 28 days are compared with 
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Figure 2. Distributions P{Bt) of the original and rescaled burstiness of invididual 
users with the same strength (left) and distributions P{ABt) of the difference in 
burstiness, defined as ABt = Bt — Bq (right). The individual users with the strengths 
Si — 200 (a), 400 (b), 800 (c), and 1600 (d) are analyzed. The numbers of users are 
correspondingly 6397, 1746, 196, and 7. 

the original distributions in the right column of Fig. [Tj Note that the possible minimum 
value of rescaled inter-event time is T/sj. We find that the rescaled inter-event time 
distributions still show the heavy tails. For the user with strength 200, the burstiness 
decreases from the original value of Bq ^ 0.202 to value Bi ^ 0.174 (weekly pattern 
removed), then dropping further to value i?28 ~ 0.104 (i.e. monthly pattern removed). 
For the most active user with strength Sj = 3197, the burstiness decays faster as T 
increases: Bq ^ 0.469, Bj ^ 0.254, and B28 ~ 0.219. However, the values of B are 
overall larger than those of the less active user. The results imply that de-seasoning 
the circadian and weekly patterns does not considerably affect the temporal burstiness 
patterns of individuals. Finally, in a limiting case of T = Tf, since ni{t) has the value 
of either or 1, all r* are the same as T/si in Eqs. ([T]) and ([s]), leading to Bxf = — 1. 

Next, we obtain the distributions P{Bt) of original and rescaled values of burstiness 
of individual users with the same strength and the distributions P{ABt) of the difference 
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Figure 3. Burstiness Bt as a function of period of T: (a) the average and the standard 
deviation of Bt obtained from the burstiness distribution in Fig. [2] (b) burstiness from 
groups with the same strength, and (c) burstiness from groups with broad ranges of 
strength. 



in burstiness, defined as ABt = Bt — Bq, see Fig. [2j The averages and the standard 
deviations of burstiness distributions for different periods of T are plotted in Fig. |3](a). 
The more active users have the larger values of burstiness, while the values of burstiness 
of the more active users decays faster (slower) than those of the less active users before 
(after) T = 7 days. The overall behavior of the distributions shows that de-seasoning the 
circadian and weekly patterns does not destroy the bursty behavior of most individual 
users irrespective of their strengths. In addition, we find some exceptional users whose 
original values of burstiness are negative, indicating for more regular behavior than the 
Poisson process, and we also find a few individual users whose values of burstiness have 
grown as a result of de-seasoning, i.e. ABt > 0. 



2.2. De-seasoning the groups of individuals with the same strength 

Here we analyze the group of individual users with the same strength, i.e. Ag = {i\si = 
s}. The averaged event rate of a group is measured by merging individual event rates, 
precisely by obtaining ^A^(t) = XlieA^ ni^t). Figure |4] shows the original and the rescaled 
event rates with T = 1 day (left) and the original and the rescaled inter-event time 
distributions with various periods of T (right) for groups with strengths s = 200, 400, 
800, and 1600. The values of burstiness decrease only slightly as T increases, but are 
smaller than those of the original burstiness, as shown in Fig. |3](b). 

The burstiness of groups of individuals with the same strength is larger than 
the average values of individual burstiness from P{B) of the same strength. For 
example, Bq ^ 0.256 for the group of strength s = 200 turns out to be larger 
than J P{Bo)dBo ^ 0.204. Regarding to this difference, we would like to note that 
the de-seasoning of individual event times by means of the averaged event rate may 
cause systematic errors due to the different circadian and weekly patterns between the 
individual and the group. For resolving this issue, the various data clustering methods, 
such as self-organizing maps, can be used to classify users' activity patterns beyond 
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Figure 4. De-seasoning of MPC patterns for groups with the same strengths: s = 200 
(a), 400 (b), 800 (c), and 1600 (d). The original and rescaled distributions of burstiness 
are plotted in Fig. [sj^b). 



their strengths and then perform the de-seasoning separately for the different groups. 

2.3. De-seasoning the groups of individuals with broad ranges of strength 

For the larger scale analysis, we consider the strength dependent grouping of users, i.e. 
groups of individual users with a broad range of strengths, denoted by Ami,m2 = < 
Si < m2}, as similarly done in [T21 |2n]- The values of ms are determined in terms of the 
ratio to the maximum strength Smax = 7911, see Table [T] for the details of the groups. 
We determine the averaged event rates of the groups and some of them are shown in 
the left column of Fig. |5j By means of the event rates, we perform the de-seasoning 
to get the rescaled inter-event time distributions, see the right column of Fig. |5j It is 
found that the values of burstiness initially decrease slightly and then stay constant at 
relatively large values as T increases, shown in Fig. [sj^c). These results again confirm 
our conclusions that de-seasoning the circadian and weekly patterns does not wipe out 
the bursty behavior of human communication patterns. 
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Table 1. Strength dependent grouping of individuals in the MPC dataset. For each 
group, the range of strength, also in terms of the ratio to the maximum strength 
Smax = 7911, the number of users, and its fraction to the whole population except for 
the user with Smax are summarized. 



group index strength range (%) the number of users (%) 
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2.4- Power spectra analysis 

In order to see clearly the effect of de-seasoning on the event rates, we compare the power 
spectra of the rescaled event rates to the original event rates. The power spectrum of 
the event rate is defined as 

2 

(5) 

where / denotes the frequency. In Fig. [6] we show the results of this comparison, 
where it is evident that with our de-seasoning methods, the circadian and weekly 
peaks of the original power spectrum are successfully removed and do not show in the 
rescaled spectra. In all the cases, ranging from the individual de-seasoning to the whole 
population de-seasoning, when T = 1 day, the circadian peak at 1 // = 1 day is removed 
while the others, i.e. the weekly and monthly peaks remain. For T = 7 days, the weekly 
peak at 1// = 7 days is removed, and so on. This behavior can be understood because 
the cyclic patterns longer than T will not be affected by de-seasoning with the period 
T. 

3. Summary 

The heavy tails and burstiness of inter-event time distributions in human communication 
activity are affected by circadian and weekly cycles as well as by correlations rooted in 
human task execution. To investigate the existence of correlations rooted in human task 
execution we devised a systematic method to de-season circadian and weekly cycles 
appearing in human activity and successfully demonstrated their removal. Here the 
circadian and weekly patterns extracted from the mobile phone call and Short Messages 
records are used to rescale the timings of events, i.e. the time is dilated or contracted 
during high or low call or SM activity of individual service users, respectively. 
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Figure 5. De-seasoning of the MPC patterns for groups of individuals with broad 
ranges of strength: groups 1 (a), 3 (b), 5 (c), 6 (d), and the whole population (e). For 
the details of groups, see Table [l] The original and rescaled distributions of burstiness 
are plotted in Fig. [sj^c). 

We have found that after de-seasoning circadian and weekly cycles driven 
inhomogeneities, the heavy tails and burstiness of inter-event time distributions still 
remain. Hence our results imply that the heavy tails and burstiness are not only 
the consequence of circadian and weekly cycles but also due to correlations in human 
task execution. In addition, we calculated the Fano and Allan factors [12] for the 
de-seasoned mobile phone communication data, which yielded further evidence of the 
remaining burstiness in both mobile phone call and Short Messages activities of human 
communication. Although beyond the scope of the current study, as the next step 
one needs to focus on describing the mechanisms of the remaining burstiness in human 
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Figure 6. Power spectra, P{f), of the original and the rescaled event rates with 
various periods of T for individual users with strengths 200 (a) and 800 (b), for groups 
with the same strengths of 200 (c) and 800 (d), and for the whole population (e). The 
circadian and weekly peaks in the original power spectra are successfully removed by 
the de-seasoning in various ways. 



communication patterns, served well by building models and analyzing their results. 
Appendix 

We also investigated the effect of circadian and weekly cycles on the heavy-tailed inter- 
event time distribution by using Short Message (SM) dataset from the same European 
operator mentioned in the main text. We have only retained links with bidirectional 
interaction, yielding = 4.2 x 10^ users, L = 8.5 x 10^ hnks, and C = 114 x 10^ events 
(SMs). We have merged some consecutive SMs sent by one user to another within 10 
seconds into one SM event because one longer message can be divided into many short 
messages due to the length limit of a single SM (160 characters) 

We apply the same de-seasoning analysis described in the main text to the SM 
dataset. Figure [7] shows the various circadian activity patterns of individual users. The 
main difference from the MFC dataset is that the activity peak is around 11 FM. This 
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Figure 7. De-seasoning SM patterns of individual users: the original and the rescaled 
event rates with period of T = 1 day (left) and the original and the rescaled inter-event 
time distributions with various periods of T (right). Individual users with the strength 
= 200 (a), 400 (b), 800 (c), 1600 (d), and 3200 (e) are analyzed. The original 
inter-event time distribution of the whole population is also plotted as a dashed curve 
for comparison. 



feature becomes evident if the averaged event rates are obtained from the same strength 
groups or from the groups with broad ranges of strength, shown in the left columns of 
Figs [To] and [TTj For the details of the strength dependent grouping, see Table |2j 

The inter-event time distributions and their values of burstiness are also compared. 
At first, the distributions cannot be described by the simple power-law form but by the 
bimodal combination of power-law and Poisson distributions as suggested by |23]. As 
shown in Fig. |9l the values of burstiness slowly decrease as the period T increases up 
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Figure 8. Distributions P{Bt) of the original and rescaled burstiness of invididual 
users with the same strength (left) and distributions P{ABt) of the difference in 
burstiness, defined as ABt — Bt — Bq (right). The individual users with the strengths 
Si — 200 (a), 400 (b), 800 (c), and 1600 (d) are analyzed. The numbers of users are 
correspondingly 1434, 344, 62, and 10. 
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Figure 9. Burstiness Bt as a function of period of T: (a) the average and the standard 
deviation of Bt obtained from the burstiness distribution in Fig. [s] (b) burstiness from 
groups with the same strength, and (c) burstiness from groups with broad ranges of 
strength. 
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Figure 10. De-seasoning of SM patterns for groups with tire same strengths: s = 200 
(a), 400 (b), 800 (c), and 1600 (d). The original and rescaled distributions of burstiness 
are plotted in Fig. [9]^b). 



Table 2. Strength dependent grouping in the SM dataset. For each group, the range 
of strength, also in terms of the ratio to the maximum strength s,„ax — 37925, the 
number of users, and its fraction to the whole population except for two most active 
users are summarized. 



group index 


strength range (%) 


the number of users (%) 





0-94 (0-0.25) 


3685116 (87.2) 


1 


94-189 (0.25-0.5) 


294181 (6.9) 


2 


189-379 (0.5-1) 


152893 (3.6) 


3 


379-758 (1-2) 


63379 (1.5) 


4 


758-1517 (2-4) 


22507 (0.53) 


5 


1517-3034 (4-8) 


7505 (0.18) 


6 


3034-6068 (8-16) 


2393 (0.057) 


7 


6068-12136 (16-32) 


316 (0.007) 


8 


12136-24272 (32-64) 


30 (0.001) 


whole 


0-24272 (0-64) 


4228320 (100) 



to 7 days, which imphes that de-seasoning the circadian and weekly patterns does not 
considerably affect the bursty behavior. 
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Figure 11. De-seasoning of tlie SM patterns for groups of individuals with broad 
ranges of strength: groups 1 (a), 3 (b), 5 (c), 7 (d), and the whole population (e). For 
the details of groups, see Table |2] The original and rescaled distributions of burstiness 
are plotted in Fig. [oj^c). 



Finally, we perform the power spectrum analysis for SM dataset in Fig. 12 and find 
again that the cyclic patterns longer than T cannot be removed by the de-seasoning with 
period T. All these results confirm our conclusion that the heavy tail and burstiness 
are not only the consequence of circadian and other longer cycle patterns but also due 
to other correlations, such as human task execution. 
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Figure 12. Power spectra, P{f), of the original and the rescaled event rates with 
various periods of T for individual users with strengths 200 (a) and 800 (b), for groups 
with the same strengths of 200 (c) and 800 (d), and for the whole population (e). The 
circadian and weekly peaks in the original power spectra are successfully removed by 
the de-seasoning in various ways. 
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